Discrete Nonlinear Breathing Modes in Carbon Nanotubes 
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We study large-amplitude oscillations of carbon nanotubes with chiralities (m,0) and (m,m) and 
predict the existence of localized nonlinear modes in the form of discrete breathers. In nanotubes 
with the index (m, 0) three types of localized modes can exist, namely longitudinal, radial, and 
twisting breathers; however only the twisting breathers, or twistons, are nonradiating nonlinear 
modes which exist in the frequency gaps of the linear spectrum. Geometry of carbon nanotubes 
with the index (m, m) allows only the existence of broad radial breathers in a narrow spectral range. 



Carbon nanotubes [H are attracting considerable at- 
tention in recent years after their discovery by lijima 0. 
They can be thought of as cylinders of carbon atoms ar- 
ranged in hexagonal grids and are thus fullerene-related 
structures. The growing interest to study carbon nan- 
otubes can be explained by their unique physical prop- 
erties and their potential for a wide range of possible 
applications. In particular, the carbon nanotubes are 
known for their superior mechanical strength Q and 
good heat conductance In addition, it is well estab- 
lished that Cgo fuUerenes can support large-amplitude 
oscillations [H] which can be excited and controlled by 
temporally shaped laser pulses 

In the continuum approximation, nonlinear dynam- 
ics of carbon nanotubes has been analyzed by several 
groups and, in particular, supersonic longitudinal com- 
pression solitons described by the effective Korteweg-de 
Vries equation have been predicted to exist in such struc- 
tures, similar to other simpler discrete lattices How- 
ever, the recent numerical modeling of more complete 
discrete model of carbon nanotubes has demonstrated [8] 
that acoustic solitons do not exist in such curved struc- 
tures, and their supersonic motion is always accompanied 
by strong radiation of phonons. 

In this Letter, we focus on the study of large- amplitude 
oscillating modes of carbon nanotubes that have the ad- 
ditional features of being nonlinear as well as discrete. 
We reveal that both nonlinearity and discreteness in- 
duce localization of anharmonic oscillations and, as a re- 
sult, the combination of both leads to the generation of 
specific spatially localized modes of the type of discrete 
breathers 0, These modes act like stable effective 
impurity modes that are dynamically generated and may 
alter dramatically many properties of carbon nanotubes. 

Discrete breathers appear in strongly nonlinear sys- 
tems, and their spatial size become comparable to the 
lattice spacing. Such discrete breathers-also called 'in- 
trinsic localized modes' or 'discrete solitons' are respon- 
sible for energy localization in the dynamics of discrete 
nonlinear lattices 
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12l |. The manipulation of the dis- 
crete breathers has been achieved in systems as diverse as 
annular arrays of coupled Josephson junctions [l3l |. opti- 



cal waveguide arrays pLJl , and antiferromagnetic spin lat- 
tices [iBl- The direct observation of highly localized, sta- 
ble, nonlinear excitations at the atomic level in the struc- 
tures such as carbon nanotubes will underscore their im- 
portance in physical phenomena at all scales. In this Let- 
ter we demonstrate that carbon nanotubes with the index 
(to, 0) support at least three types of discrete breathers, 
whereas nanotubes with the index (m, 0) support only 
one type of discrete breather. 

The structure of a carbon nanotube is shown schemat- 
ically in Fig. 1. In static, the nanotube is characterized 
by its radius R and two step parameters hi and In 
each layer, the nanotube has to atoms separated by the 
angular distance Acf) = 2n/m, so that hi and h2 define 
alternating longitudinal distances between the transverse 
layers. We consider such dynamics of the nanotube that 
all atoms in one transverse layer have identical displace- 
ments. In this case, the carbon nanotube can be modeled 
by an effective one-dimensional diatomic chain, where 
the coordinates of atoms {pn.i cos ipn,i, Pn,isiaipnj, Zn) 
are defined by the equations: pn.i = i? + r„(t); ipn,i = 

A(/)(Z - 1) + Mt), Zn,l = (2/c - 2)(/li + h2) + Unit) 

for 71 = 4fc - 1; ipnj = A(f){l - 1) + A0/2 + 0„(i), 
Zn,i = {2k - 2){hi + /12) + hi + Un{t) for n = Ak; Lpn.i = 

A0(/-l)-HA0/2-|-(/.„(t), Zn,l = {2k-l){hi + h2)+Un{t) 

for n = 4fc -t- 1; and Lpn.i = A0(; - 1) + Zn.i = 

{2k - l){hi + + /ii + Un{t) for n ^ Ak + 2, where 
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Figure 1: Schematic of a carbon nanotube with the index 
(m, 0). Below: an effective one-dimensional diatomic chain. 



2 



the index n = Ak + i, {k = 0, ±1, ±2, i = -1, 0, 1, 2) 
stands for the number of the transverse atomic layer, the 
index I — I, ...,m marks an atom in the transverse layer, 
r„ (t) is a relative change of the radius of the n-th trans- 
verse layer, is the angle of rotation of the atoms in 
the layer, and u„ is a relative longitudinal displacement 
of the atoms from their equilibrium position. In the static 
case, r„ = 0, = 0, and Un = 0. 

In this case, complex three-dimensional dynamics of a 
nanotube can be reduced to the analysis of an effective 
one-dimensional model with two atoms per its unit cell, 
and Hamiltonian of this model can be written in the form, 

+^2n+i+^(y2n-2; y2n-i; y2n; y2n+l)], (1) 

n 

where = M[rl + {R + rn)^(pl + ul]/2 is the kinetic 
energy , M is the mass of a carbon atom (the total energy 
of the nanotube is mH), and the vector y„ = (r^, (/)„, w„). 
The interatomic potential can be written in the form, 

^(yi,y2,y3,y4) = V^(xi,X2) + F(X2,X3) + V^(X2,X4) 
+ [/(x2,X3,X5) + [/(X6,X3,X2) + t/(x6,X3,X5) 
-|-[/(xi,X2,X3) + C/(xi,X2,X4) + t/(x3,X2,X4) 
+W^(X6 , X3 , X2 , X5 )-|-PF(x6 , X3 , X5 , X2 )-|- VK(X2 , X3 , X6 , X5 ) 
+ W^(X1, X2, X4, X3)+M^(X1, X2, X3, X4)+VF(X3, X2, Xl, X4), 

where the coordinate vectors x^ = {xi,yi,Zi), i — 
1,2,.. .,6 are defined as: xq — {R + ri)cos(A0 + 
ye = (i? 4- ri) sin(A0 + 0i), zg = -/12 + mi, xi = 
(i? + r2) cos((/)2), yi = (i? + r2)sin((/)2), zi = U2, 
X3 = (i? -I- r2)cos(A<?!) + 02), ys = [R + sm{A(j3 + 

(1)2), 2:3 = M2, X2 = (i? + rs) C0s(A(^/2 -I- (/)3), 2/2 = 

(R + r3)sin(A</./2 + ^3), Z2 = hi + U3, X5 = {R + 
r3)cos(3A0/2 -I- 03), ys = (i? + r3) sin(3A0/2 + ^3), 
Z5 = hi + U3, X4, ^ {R + r4)cos(A0/2 + (f)^), 7/4 
(-R + r4) sin(A(/)/2-|-(}()4), Z4 = hi + h2 + U4. The potential 
F(xi,X2) = D{exp(-a[p - po]) - 1}^, p = |x2 - Xi|, 
describes a change of the deformation energy due to in- 
teraction between two atoms with the coordinates Xi 
and X2. Potential C/(xi,X2,X3) — (cos(/3 + 1/2)2, 
where cos(p — (vi, V2)/(|vi| • |v2|), and vi = X2 — xi, 
V2 = X3 — X2, describes the deformation energy of the 
angle between the links X1X2 and X2X3. Finally, the po- 
tential Ty(xi,X2,X3,X4) = et(l — cos(/)), where cos0 = 
(ui,U2)/(|ui| • |U2|) and Ui = (x2 - Xi) x (xa - X2), 
U2 = (x3 — X2) X (x4 — X3), describes the deforma- 
tion energy associated with a change of the effective an- 
gle between the planes X1X2X3 and X2X3X4. We take 
the mass of carbon atom as M = 12mp, where nip if 
the proton mass, the length po = 1.418A, and energy 
D = 4.9632 eV. Other model parameters such as a, e^, 
and et can be determined from the phonon frequency 
spectrum of a plane of carbon atoms. 

A flat plane of the carbon atoms (graphene) is a special 
case of carbon nanotubes in the limit i? — > 00, when 



hi = Po/2 and /12 = Po- For such a plane the motion 
equation splits into the equations for longitudinal and 
transverse motion. The corresponding Hamiltonian takes 
a simpler form, 

H = T.^\m{uI^_i+uIJ + Vi{p2n) + V2{p2n-l)], (2) 
n 

where p„ = m„+i — u„, the potentials Vi{w) = 
D[exp{-aw) - 1]^, T/2H = 2D{cxp{-a[a+{w)'^/^ - 
po]) ~ ]V + 2£.jQ_(w)/a+H + l/2]2 + Ae^iw + 
po/2)/ v^cM~(w) — 1/2]^, and function a±{w) = {w + 
Po/2)2±3p§/4. 

The stiffness parameters of the potentials are Ki = 
V('{0) = 2Da^ and K2 = V^'{0) = Da^ + 27ej2p2. 
After linearizing the equation of motion following from 
the Hamiltonian we obtain the dispersion relation 
for longitudinal phonons, uj±{q) — {{Ki + K2 ± [{Ki + 
K2f - 2KiK2{l - cos2g)]i/2)/M}i/^ which are de- 
picted in Fig. [2] The frequency spectrum consists of the 
acoustic [w_(0), w_(7r/2)] and optical [ lli+{-k/2), uo+{Q) \ 
bands where w_(0) = 0, cj+(0) = ^2{Ki+ K2)/M, 
lj-{-k/2) = sj2K2lM, and w+(7r/2) = ^2Ki/M. 

The edge frequencies of the optical band can be 
estimated from the experimental data: uj+{'k/2) k, 
1200 cm-i and u;+(0) ^ 1600 cm"! [l3|. These val- 
ues allow us to determine the stiffness parameters, Ki — 
508.98 N/m and K2 = 395.87 N/m, and find the 
maximum frequency of acoustic phonons, u)^{ti/2) = 
1058.3 cm^^. Knowing the values of Ki and ^2, we 
then find other parameters, a = 1.7889 and = 
1.3143 eV. The value of the torsion potential et can be 
evaluated from the maximum frequency of the transverse 
oscillations of a plane carbon lattice. For et — 0.2 eV, we 
find the value 570 cm~^. In order to find the parameters 
R, hi, and /12, we should solve the minimum problem 
Z(0, 0, 0, 0) — > mini^^/ij_ft,2 . The resulting value of energy 
is then used as the minimum value. 

A simple form of the Hamiltonian ^ allows to ob- 
tain analytical results for the nonlinear dynamics similar 
to the case of diatomic lattices. These results allow to 
predict the existence of discrete breathers with the fre- 
quencies below the lowest optical frequency of the longi- 
tudinal phonons [see Fig.[2l^c)]. The form of this breather 
is shown in Fig. [3l The breather is characterized by the 
frequency w, energy E, width of the localization region, 
L, and the chain extension A. The breather frequency is 
inside the band [1162, 1200] cm^^ near the lowest edge 
of the longitudinal optical oscillations. For decreasing 
Lo, both E and A grow monotonically, and the breather 
width decreases. 

Hamiltonian ([T]) defines the motion equations of the 
system. We have studied these equations numerically and 
revealed that they support the existence of three types 
of strongly localized nonlinear modes- discrete breathers. 
The first type, longitudinal breathers, also exists in pla- 
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Figure 2: Spectral density of thermal oscillations of a carbon 
nanotube (10,0) for the temperature (a) T = 300K and (b) 
T = 30K. (c) Dispersion curves of the phonons: acoustic (1) 
and optical (2) longitudinal phonons, acoustic (3) and optical 
(4) rotation phonons, and acoustic (5) and optical (6) radial 
phonons. For comparison, red dashed lines show the disper- 
sion curves of longitudinal oscillations of a plane of carbon. 
Green color marks the frequency spectrum of breathers. 



nar carbon structures, such breathers exist in the fre- 
quency range [1162, 1200] cm^^. The second type, ra- 
dial breathers, describes transverse localized nonlinear 
modes with the frequency band [562, 580] cm~^. The 
third type, twisting breathers, characterizes localization 
of the torsion oscillations of the nanotube with the fre- 
quencies [1310, 1477] cm"-^. The frequency spectra of the 
breathers are shown in Fig. [51 

For a flat plane of carbon atoms, the longitudinal 
breathers are nonlinear modes [see Figs.[3]^a,c)], and they 
are exact solutions of the nonlinear motion equations. 
However, in the case of a curved geometry, the longitudi- 
nal breathers become coupled to transverse linear modes, 
and they always emit some radiation. This radiation is 
defined by the curvature of the nanotube and its index 
m. Therefore, the longitudinal breathers are not genuine 
nonlinear modes of carbon nanotubes, and they possess a 
finite lifetime which however may exceed a few hundred 
of picoseconds. 

The second type of discrete breathers we found is as- 
sociated with the localization of transverse radial oscil- 
lations of a nanotube. Example of this radial breather 
in the nanotube (10,0) is shown in Figs. [3^b,d). Local- 
ized out-phase transverse oscillations of the neighboring 
atoms lead to localized contraction and extension of the 




Figure 3: Example of a localized nonlinear mode of the lon- 
gitudinal oscillations described by the Hamiltonian ^ (a), 
(b) (frequency u> = 1164cm~^, energy E — 5.71eV, width 
L = 16.3, the chain extension A — 0.45A) and example of a 
radial breather (c), (d) describing localized transverse oscilla- 
tions of a nanotube (10,0) with the frequency lo = 579.6 cm~^, 
energy E = 0.99853eV, and width L = 32.7. Shown are (a) 
and (c) the averaged (in time) energy distribution En in the 
chain, (b) the atom displacements Un, and (d) transverse dis- 
placements r„. In sections (b) and (d) blue lines show the val- 
ues averaged over the period, red - maximal displacements. 
Radiation of longitudinal waves by a radial breather is clearly 
visible in section (d). 



nanotube. Such transverse oscillations become coupled 
to the longitudinal oscillations and, therefore, the radial 
breathers radiate longitudinal phonons. As a result, the 
radial breathers are also not genuine nonlinear localized 
modes of the carbon nanotubes, and they decay slow by 
emitting small-amplitude phonons. The lifetime of these 
breathers can be of the order of several nanoseconds. 

The third type of localized mode is a twisting breather, 
or twiston, associated with the torsion oscillations of the 
nanotube. In a sharp contrast to other two breathing 
modes, the twisting breather is an exact solution of the 
motion equations of the nanotube, and it does not radiate 
phonons. An example of this genuine discrete breather 
is shown in Fig. 3) In the localized region of this mode, 
the nanotube is expanded transversally being contracted 
longitudinally. The twiston has a broad frequency spec- 
trum, and its energy, amplitude of the transverse exten- 
sion (see Fig. [5]), and the amplitude of torsion oscillations 
all grow with the frequency. The breather width changes 
monotonically, and for the frequencies ui < 1450 cm^^ it 
becomes comparable with the lattice spacing, so that the 
breather becomes a highly localized mode. 

Next, we analyze thermal oscillations of a carbon nan- 
otube by employing the Langevin equations. We find 
that for low temperatures (T = 30K) the oscillations 
are mostly linear and the frequency density does not dif- 
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Figure 4: Dependence of (a) transverse r„, (b) torsion (j)„, 
and (c) longitudinal it„ displacements of the twisting discrete 
breather (frequency uj — 1367 cm~^, energy E — 1.6465 eV, 
width L — 1.55) on the time t. 
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Figure 5: Dependence of total energy E (curve 1) and width L 
(curve 2) of the twisting discrete breather on the frequency uj. 
Vertical line show the edge of the linear spectrum of optical 
torsion oscillations in the (10,0) nanotube. 

fer much from the density of linear modes, as demon- 
strated in Fig. El^b). However, for higher temperatures 
(T = 300K) the spectral density acquires specific features 
associated with the generation of nonlinear modes, as can 
be seen in Fig.[2l^a) where we observe thermal oscillations 
with the frequencies inside the linear spectral gaps which 
can be associated with discrete breathers. Indeed, the 



larger contribution of these nonlinear modes is for the 
torsion oscillations of the twisting breathers, which are 
stable and have the largest frequency spectrum. 

We have carried out the similar nonlinear analysis for 
the nanotubes with the index (m, m) and revealed that 
in this nanotube can support only one type of breathers, 
a radial breather with a very narrow frequency spectrum 
[430.5, 436] cm~^ near the upper edge of the frequency 
band of radial phonons. However, this radial breather is 
not an exact solution of the nonlinear motion equations, 
and radiation of small-amplitude linear waves leads to a 
decay of the breather. As a result, the existence of nonlin- 
ear localized modes depends crucially on chirality of the 
carbon nanotube, so that genuine discrete breathers are 
expected to exist in the nanotube with the index (m, 0). 
Existence of twisting breathers is due to the anharmonic 
interaction potential, and the large spectral gap in the 
frequency spectrum of torsion phonons. However, the 
radial long-lived nonlinear modes can appear in the nan- 
otubes with any type of chirality. 

In conclusion, we have revealed that carbon nanotubes 
can support spatially localized large-amplitude stable 
nonlinear modes in the form of discrete breathers, and we 
have analyzed the existence and stability of three types 
of breathers. A novel type of such highly localized dis- 
crete modes -twistons- is associated with the energy self- 
trapping of torsion oscillations of the carbon nanotubes. 
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